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ABSTRACT 

We investigated the intracellular delay effect on the 
stability of the endemically infected steady state by 
analyzing a nonlinear ordinary differential equation 
model of hepatitis B virus (HBV) infection that 
considers the interaction betweena replicating virus, 
hepatocytes and the cytotoxic T lymphocytes (CTL). 
We gave a criterion to ensure that the infected steady 
state is asymptotically stable for all delays. A critical 
delay below which the CTL (immune control 
mechanism) can be significantly helpful in controlling 
the HBV infection even when the basic reproduction 
number is high is allowed in the analysis. 

KEYWORD: Delay differential Equation, Critical 
Delay, Hepatitis B Virus, immune system 

I. INTRODUCTION 

Hepatitis B virus (HBV) infection is a major health 
problem worldwide and despite the presence of 
hepatitis B vaccine, new HBV infections remain 
co mm on. HBV is a hepatotropic double stranded 
DNA virus belonging to the hepadnaviridae family. 
One of the five hepatitis viruses, HBV causes acute 
and chronic hepatitis in humans. Acute HBV infection 
causes only mild symptoms. Majority of infected 
adults successfully clear the virus and acquires 
lifelong immunity. It takes about 1 to 6 months from 
the time of infection for the disease to manifest itself 
and early symptoms include nausea and vomiting, loss 
of appetite, fatigue, with dark urine and light stools 
follows. Only about one percent of patients infected 
with hepatitis B virus die due to liver damage in their 
early stage (Poynard, 2002). The risk of becoming 
chronically infected depends on the age at the time of 
infection. More than 90 percent of new bom, 50 
percent of children, 5 percent of adults infected with 


HBV develop chronic hepatitis (WHO, 2014). Those 
who are unable to produce an effective immune 
response allow the vims to replicate for long periods 
in their liver causing chronic hepatitis B vims, (HBV) 
infection, cirrhosis of liver and hepatocellular 
carcinoma (HCC) (Yen, 2002). Chronic HBV 
infection is shown to be 0.8% in pregnant women 
with the risk of vertical (mother to child) transmission 
(Trehanpati et al, 2013). In neonates, HBV infection 
has a linkage with strong presence of Tregs (T 
regulatory cells) and defective CD8 + T cells pool to 
produce interferon QFN) — y. Associated with 
decreased CD8 + T cell dysfunction is T cell receptor 
(TCRQ chain. Due to persistent intrauterine exposure 
of the viral antigens early in embryonic development 
there could be TCR £ leading to immune tolerance to 
HBV antigens in the newborns positive for hepatitis B 
surface antigen (HBsAg +ve). Finally T cell tolerance 
to HBV antigen due to HBV infection may probably 
develop to chronicity in the newborn. To restore 
acquired immunity and better production of HBV 
specific antibodies, vaccination may be necessary 
(Trehanpati et al, 2013). 

Two disease states, not necessarily static, are possible 
after seroclearance of HBeAg and they are inactive 
carrier (IC) state and chronic hepatitis B (CHB) state 
(Pita et al, 2014). The duration of clearance (reactive) 
phase last from months to years. In the inactive carrier 
(IC) state, seroconversion of hepatitis B e antigen 
(HBeAg) to hepatitis B e antibody (HBeAb) occurs, 
HBV DNA becomes non detectable or at low level 
and ALT is usually normal, reflecting very low or no 
replication of HBV and mild or no hepatic injury 
(European association for the study of liver (EASL), 
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2012). The inactive carrier state may last for years or 
even lifetime. Patients in this state can have 
spontaneous resolution of hepatitis B and develop 
HBsAg, but a portion of them may undergo 
spontaneous or immunosuppressant - induced 
reactivation of chronic hepatitis B. This bring in the 
second state which is chronic hepatitis B and is 
usually associated with elevated ALT, high level of 
DNA, moderate to severe liver histological activity, 
and with or without HBeAg seroconversion (Pan and 
Zhang, 2005). Differentiating between inactive carrier 
(IC) and chronic hepatitis B (CHB) status is very 
important in clinical practice, because it has an 
implication in the follow-up, management and 
prognosis (EASL, 2012). 

The HBV virion genome is circular and 
approximately 3.2kb in size, consisting of DNA that is 
mostly double stranded. It has compact organization, 
with four overlapping reading frames running in one 
direction and no non coding regions. The minus 
strand is unit length and has a protein covalently 
attached to the 5' end. Thus, neither DNA strand, is 
closed and circularity is maintained by cohesive ends 
(Strauss, 2002). The four overlapping open reading 
frames (ORFs) in the genome are responsible for the 
transcription and expression of seven different 
hepatitis B proteins. The translation of these proteins 
is through the use of multiple in-frame start codons. 
The HBV genome also contains parts that regulate 
transcription, determine the site of polyadenylation 
and a specific nucleocapsid. The genomic 
arrangement of the hepatitis B virus family makes it 
unique among viruses. The unusually packaged may 
indicate that the method of replication employed by 
HBV is not of conserved DNA replication (Graces, 
HBVP, 2003). 

Since HBV is a noncytopathic virus (Ilan, 2002), that 
is cells infected by HBV will not be killed by virus 
directly, cellular function and life span of HBV - 
infected hepatocytes are almost the same as that of 
the uninfected cells in vitro (Kangxian, 2006). The 
death rate of noncytopathic virus infected cells in the 
absence of immunity equals that of uninfected target 
cells (Woderz, 2005). The lifespan of HBV infected 
cells varies greatly in vivo which is mainly due to the 
strength of the anti HBV CTL response (Nowak and 
Robert, 2000). CTL will not only kill but cure the 
infected hepatocytes by nonlytic effector mechanism 
(Bartolett a and Ferrari, 2003; and Guidotti, 2003). In 
a study by Yang et al (2009) it was observed that, 


induction of a CD 8 + T cell response depends on the 
presence of CDA + T cells. That each of these 
effectors, CDA + and CD8 + T cells, Nk cells, fas, IFN 

- gamma (IFN - y), IFN - alpha/beta receptor (IFN - 
a/P - R), and IFN receptor 1 (INF Rl), was required 
to eliminate the transcriptional template from the 
liver. These results are consistent with a model in 
which CD4 + T cells serve as master reguletors of the 
adaptive immune response to HBV; CD 8 + T cells are 
the key cellular effectors mediating HBV clearance 
from the liver, apparently by a fas - dependent 
perforin independent process in which NK cells, IFN 

- y, INF Rl and IFN - a/p - R play supportive roles 
(Yang et al, 2009). Immediately after infection, 
innate immunity limits the spread of the pathogen and 
initiate efficient development of an adaptive i mm une 
response. Innate host responses during the early phase 
of viral infection are mainly characterized by the 
production of type 1 interferon (IFN) - a/p cytokines 
and the activation of natural killer (NK) cells. 
Production of type 1 IFNs can be triggered directly by 
virus replication through cellular mechanisms that 
detect the presence of viral RNA or DNA 
(Alexopoulou et al, 2001; Lund et al, 2003; Heil et al, 
2004), while NK cells are activated by the recognition 
of stress - induced molecules and /or the modulation 
of the quantity of major histocompatibility complex 
(MHC) - class 1 molecules on the surface of infected 
cells (Moretta et al, 2005). 

II. RELATED WORKS 

Wu and Su (2014) considered a model with 
intracellular delay and nonlinear infection rate of 
saturated functional response. By constructing 
suitable Lyapunov functions and using LaSalle 
invariance principle, the proved that the infection free 
equilibrium is globally asymptotiacally stable when 
R 0 < 1. And if R 0 > 1, there exists an infection 
equilibrium and they get its global asymptotic 
stability by constructing Lyapunov functional. Wang 
and Tian (2013) developed and investigated a model 
that studied the global stability for a delayed HBV 
infection with CTL immune response. They showed 
that the global dynamics is determined by the basic 
reproduction number R 0 and CTL immune response 
reproduction number R t . 

Huang et al (2010) studied the global properties of a 
class of delay differential equations model for virus 
infection with nonlinear transmissions. Their model 
gives a more complex and general infection process, a 
general nonlinear contact rate between target cells and 
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viruses and the removal rate of infected cells are 
considered, and two constant delays are incorporated 
into their model. These incorporated delays describes 
(i) the time needed for a newly infected cells to start 
producing viruses and (ii) the time needed for a newly 
produced virus to become infectious (mature). Wang 
et al (2007) developed and analyzed the HBV 
infection in a diffusion model confined to a finit 
domain, induced by intracellular time delay between 
infection of a cell and production of a new viral 
particle. Rezounenko (2017) investigated a class of 
reaction diffusion virus dynamics model with 
intracellular state dependent delay and a nonlinear 
infection rate functional response. 

III. Formulation of the Model Equations 
3.1 The Existing Model 

We start our model formulation by introducing the 
model by Chinebu et al (2017) and finally incorporate 
the delay parameter to obtain the delay model. 

dx Pa 3 xy 

— =xp- ai x -+ b 4 yz 

dt /i 

dy /?a 3 xy 

~dt = ^ y + —[i - aiy ~ ^ bl + b ^ yZ ( ' 3 ' 1 ) 

dz 

— = y + a 4 yz - b 4 z. 

Equation (3.1) contains three variables, that is, 
uninfected hepatocytes (x), infected hepatocytes (y), 
and CTL response (z). Uninfected, susceptible 
hepatocyte are assumed to be produced at a constant 
level, ip, the natural death rate of both infected and 
uninfected hepatocytes are assumed to be the constant 

(3a 3 xy _ s assume( j t0 b e q le l eve l a t which 

uninfected hepatocytes are being infected due to 
interaction with infected hepatocytes and <y is the 
level at which infected hepatocyte proliferate. The 
elimination of virus is through the response of the 
cytotoxic T lymphocytes (CTL) which can be 
activated in two different pathways, either by killing 
the infected hepatocytes or by clearing the virus from 
within the hepatocyte without killing it. Infected 
hepatocytes are assumed to be killed by the CTL 
response at a rate b 2 yz and be cleared from within by 
the CTL response at a rate b 4 yz. a 3 y is the 
production level of free virus from infected 
hepatocytes and the free virus particles are removed at 
a rate /.(.Two terms, y and a 4 yz were used to describe 


CTL rapid growth and reproduction of new cells 
(proliferation), where y represents antigen 
independent proliferation (primary immune response) 
anda 4 yz represents antigen dependent proliferation 
(secondary immune response), y is small, and is the 
response to the first exposure to antigen which is 
mediated by lymphocytes, called naive lymphocytes 
that are seeing antigen for the first time and they are 
“immunologically inexperienced”. We also assume 
that CTL decay at constant rate b 4 (Chinebu et al, 
2017). 

3.2 Delay Differential Equation 

A delay differential equation is a type of differential 
equation whereby the derivative of the unknown 
function at a certain time is given in terms of the 
values of the function at previous time and is also 
known as time delay system. Notice that we did not 
incorporate any time delay in HBV proliferation in the 
model presented in equation 3.1. Between the 
infection of a cell and the viral particles emission, 
virus production may lag by an intracellular time 
delay r. The delay or lag represents incubation period, 
thereby accounting only for the time required for the 
production of new virus particles. However, Rotich 
and Lag at (2014), shows that when the drug efficacy 
is less than perfect, the intracellular delay will be of 
crucial importance because the delay depends on the 
live span of the liver cell, that is to say that the 
maximum length of time that the cell will survive 
after infection is considered. The incorporation of 
intracellular delay is assumed to give a clearer view in 
the interpretation of the analytical result and this has 
important implications on therapeutic options and 
drug development. 

Gourley et al (2008), developed a delay differential 
equations on HBV infection which explicitly 
incorporate a time delay in virus production. For their 
model, the existence and the component values of the 
endemic steady state are explicitly dependent on the 
time delay. They further stated that exposed cells 
begin shedding virions after r unit of time, 
representing the time required to construct, transcribe 
and translate the viral genome. The episomal viral 
genome, construct and then release mature virions. 
There work extended the ability of models to describe 
the relevant biological process and addressed the 
implication of ignoring the intracellular delays that 
are part of the viral life cycle. However, they did not 
take into cognizance that cytotoxic T lymphocytes 
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play an important role in the immunological 
mechanism as will be seen in our later discussion. 
Consider a system that can be in either one of two 
states Q (infection of a cell) and R (viral particle 
emission). Suppose that at timet = 0, the system is in 
state Q. suppose that when event p takes place, the 
system switches to state R. Modifying the model 
given by (3.1), we have a new model which assumes a 
delay r, from the time of infection until the production 
of new virus particles. We define a translation 
operator in order to describe the delay in 
mathematical terms. 


Where functions F lt F 2 and F 3 depends on x, y, z 
and their delays. Therefore, F v F 2 and F 3 can be 
identified as follows: 


Ft(x, y, z, x 1 ,y 1 ,z 1 ) 


= i p — a 4 x — 


/?q 3 *y 

F 


+ b 4 yz 


F 2 (x,y, z, x 1 ,y 1 ,z 1 ) 

pa 3 x 1 y 1 

= O + — -% y 

[l 


- Oi + b 2 )yz 


(3.4) 


F 3 (x,y,z,x 1 ,y 1 ,z 1 ) = y + a 4 yz - b 4 z. 


For a real - value function p and for r > 0, we define 
the translation operator T T by 
(T T p)(t) = p(t — t), with time delay t > 0 and p 
E R n . 


Since the virus must have contact with the liver cells 
before they will be infected, we let 

*i = T t x, y 4 = T T y, z 4 = T t z 

for notional convenience and this gives the new 

system with delay as 

dx fia 3 xy 

— = ill — a.x -1- hiyz 

dt p 


dy Ba 3 Xiyi 

— = +- ayy - (b 4 + b 2 )yz 

dt fi 


(3.2) 


dz 

— = y + a 4 yz - b 4 z. 


^ a 3*i y i , indicates a finite time lag between infection 

of a hepatocyte and production of new virus particles 
in the equation that determine the rate of y and this 
describes the dynamical variable y at time t that 
depends on the term evaluated at earlier time t — r. 
Since there are three, possible steady state in the 
previous section and these states are too complex to 
analyze, we therefore, discuss only the linear ability 
of the most concerned state - uninfected state. System 
(3.2) is a special case of the general system: 
dx 

— = F 4 (x, y, z, x 4 ,y x ,Z]) 


— = F 2 (x, y, z, x lf y lt z 4 ) (3.3)' 


dz 

— = F 3 (x,y,z,x 1 ,y 1 ,z 1 


) 


We then use define the steady state for the delay 
system having that a point (x, y, z) is called a steady 
state of the system (3.3) if it is a constant solution of 
the equations. 

F 4 (x, y, z, x 4 , y 1( z 4 ) 

F 2 (x,y,z,x 1 ,y 1 ,z 1 ) (3.5) 

F 3 (x, y, z, x 1 ,y 1 ,z 1 ) 

Propositionl: If a ll a 3 ,b ll b 2 ,b 4 ,\p,^ l Y,f3,p> Q, 
then the steady state of system (3.2) is the same as the 
steady state of system (3.1). 

Proof: From equation (3.4), we have that 

4 0 K,y,z,x 1 ,y 1 ,z 1 ) = /i(x, y, z) 

F / A A A A A AN r / A A AN 

2 (x,y, z, x 1 ,y 1 ,z 1 ) = f 2 (x,y,z) 

3 (x, y, z, x 1( y lf zj = f 3 (x, y, z) 

This simply shows that the steady states of system 
(3.2) and (3.1) are the same. Considering a 
neighborhood close to the steady state solution and 
letting 

x = x + X, y = y + Y, z = z + Z, 
x 4 = x + X x , y 4 = y + y±, z 4 = z + Z 4 


after which we expand F v F 2 andF 3 in the Taylor 
series about the point (x,y,z) and retain only the 
linear terms, to obtain 


dX 

dt 


- a ± 


Pas9 


F 


X + 


/3a s x 

p 


+ b 4 z 


Y + [b 4 y]Z 


dY 

dt 


Pa 3 y 


F 


Zi + 


pa 3 x' 
- F - 


n 


+ [<■-«!- (h 4 + b 2 )z]Y 
~ [(^1 + b 2 )y]Z 


(3.6) 


dZ 

— = [a 4 z]Y + [a 4 y - b 4 ]Z 
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To linearize (3.6) we let 
r dx- 


dt 

dY 


dt 

dZ 


dt 


x(t) 

y(0 

z(t) 


Let 


-X- 


-X- 

y 

= H and 

y 

z- 


z- 


(0 = H(t) 


ff-f 


Adt 


In H = At + c 


We take exponential of both sides to get 

^ln H _ gAt+c 

=> H(t) = e At . e c 


But, in biological population (Malthusian growth 
model), the major underlying assumption is that the 
rate of change of population is proportional to the 
population itself. 

Thus we let 

dH 

—— = AH. 
dt 

where the parameter A is known as the natural rate of 
increase and is a positive constant. We now solve the 
differential equation using the method of separation of 
variables to obtain 

dH 

— Adt. 

H 


If we let e c = H 0 we have 

Hit) = H 0 e Xt . 

This predicts exponential population growth as t 
increases. However, due to the time lag between the 
infection of a hepatocyte and viral particle emission 
from the hepatocyte, it may be more realistic to 
assume that the instantaneous rate of change of 
population growth is actually dependent upon the 
population at some fixed amount of time r in the past 
(Cain and Raynolds, 2006). Observe that the rate of 
change of H at time t is affected by the value of H at 
time t — r. Therefore we have 

H(t) = H 0 e A(t_T) 

=> Hit) = H 0 e Xt .e~ Xr 

At t = 0 which is assumed to be the time of infection 
of a hepatocyte, we have 

Hit) = H 0 e~ Xr 


xit) 

yit) 

z(t) 


= H n e 


-At 


where H 0 is a constant vector (Mac Donald, 1970) and 
we have 


-Xi 


-x- 

Yi 

= e Xt 

y 

Zl. 


z- 


which transforms system (4.22) into a linear system as stated below 


vdXi 

r 

dt 


dY 


dt 


dZ 


Ldt-I 

- 


d 

/?q 3 ye~ At 

H 

0 


Then, we arrive at the matrix 

/?<* 3 y\ 


( P a zy\ (Pa 3 x \ 

-l a ‘ + — ) -{—~ Wz ) 


f^a 3 x 


b±y 


< +-:-% - ib ± + b 2 )z -(^ + b 2 )y 


H = 


-( 


CL-t + 

d 

d 

0 




n 


(Pa 3 x 


a 4 z 


a 4 y — b 4 _ 


( + 


V n 

fda 3 xe~ Xt 


-m) 


bi y 


d 


-a 1 - ib t + b 2 )z ~ib x + b 2 )y 


X 

Y 

LZJ 


(3.7) 


a 4 z a 4 y — h 4 _ 

Equation (3.7) is the coefficient matrix of the linear system having A as the eigenvalue of the matrix H. At this 
point, we investigate the stability of the delay system in terms of the eigenvalues of the matrix H. 
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Delay differential system requires information on the entire interval [t 0 — r, t 0 ]. To know the rate of change at 
t 0 , we need x(t 0 ),y(t 0 ),z(t 0 ),x(t 0 — r),y(t 0 — t) and z(t 0 — t). According to Burton (1985), it is well 
known that the solution to s'(t) = s(t) 2 diverge to infinity in finite time. Solutions to the delay differential 
equation s'(t) = s(t — r) 2 , however, are continuable for all time if r is positive for t. Also,Halanay (1966), 
states that if ( x*,y*,z*) is the steady state solution of system (3.3), we say that (x*,y*, z*) is stable if for every 
a > 0 there is a S > 0 such that, for every solution (x, y, z), if 

[x(t 0 ) - X *] 2 + [y(t 0 ) - y*] 2 + [z(t 0 ) - z*] 2 < d 2 (3.8) 

For some toe^ — r, tj, then 

[x(t 0 ) - X *] 2 + [y(t 0 ) - y*] 2 + [z(t 0 ) - z *] 2 <£ 2 v t>t ± . 


If, (x*, y*, z*) is not stable, we say that it is unstable. If (x*, y*, z *) is stable and if d can be chosen so that (3.8) 
implies that (x(t),y(t),z(t)) -> ( x*,y*,z *) as t -> oo, we say that (x*,y*, z*) is asymptotically stable. From 
Halanay (1966), we remark that the steady state (x,y,z) is stable if no eigenvalue of H has positive real part. 
We compute the characteristic equation by seeking for exponential solution of the delay differential system to 
arrive at a transcendental equation of the form 

dj (A) + kj(X)e~ Xx = 0, 

Where r is the length of discrete delay added, d and k are polynomial in A and j = 1,2,3. 


Proposition!: Suppose 


a ± , a 3 , b x , b 2 , b 4 , y, X /?, ii, ip > 0 


1. If a x + 


(& 1 +& 2 ):k 


stable for r > 0 . 
Cbi+b 2 )y 


2. If a x + 


b 4 


> <( + ^ a3 ^ , then the steady state (^/ a , 0 , ^/ h ] of the system (18) is asymptomatically 
a lM \ 1 1 I 

< ( + then the steady state (^’/a 1 > 0 > ^/b^J ^ ie s Y stem (3-2) is unstable for r > 0 . 


Proof: 

At steady state, the characteristic equation of the delay differential equation will have the form 
= A 3 H~ (di "I - k^e ^ T> )A 2 + (d 2 + /c 2 £ ^ T> )A "I - (^3 "I - k%e = 0 (3-9) 


which can be written as 

t/ 2 (>l, t) = A 2 + d^A 2 + d 2 A + + (M 2 A- /c 2 A + k^e T =0 


(3.10) 


( a. pa * y \ i 

pa^ye^ 

H 

0 


//?a 3 x 


<T + 


V n 

(3a 3 xe~ At 


m) 


bi9 




- a 1 - (hi + b 2 )z - A ~{b 1 + b 2 )y 


CL4.Z 


a 4 y — b 4 — A 


= 0. 


Since the steady state about which we have linearized is stable in the absence of delay (Chinebu et al, 2017), 
that is when r = 0 ,which simply implies that there is no delay and that all the roots of the polynomial have non 
positive real parts, therefore, proposition 2 is proved as in (Chinebu et al 2017). We further consider the 
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situation where r > 0 and therefore substitute the steady state (^ > /a 1 > 0< ' nt0 system (3.7), after which 

we obtain the characteristic equation as; 

—a ± — A 


pcwp | b x y 

CL-^H b 4 


0 

0 


„ , P a 3^ e At n , , % y i 

<T H- a i ~ (, b i + b 2 ) T— A 

b 4 

a 4 y 


—b 4 — A 


= 0 


a+«n (a+ ai + f a+bt) = 0 . 

V & 4 ail* J 

Therefore, matrix H has three eigenvalues, at steady state /a 4 < 0< Vb ^and they are A 1 — —a 1 < 0 ,A 2 = 

—b 4 < 0, and A 3 satisfying the equation 

1 , , (&i + b 2)Y 7 Pa 3 i/je~ rA 

b 4 a 1 /j. 

At this point, we introduce a function v(t ) to enable us solve for the position of A 3 

^ , , C b i + b 2)Y 7 Pa 3 x/je- rt 

v(t) — t + Ci'i H-:-f-, t G M 

b 4 a x ii 

and this enable us to find the location of the eigenvalue A 3 . By differentiating v(t) with respect to t, yields 

, , P^xe~ xt 

a^u 

which is always positive, we also observe that 

lim v(t ) = —oo, lim v(t ) = oo. 

t —> — OO t >00 

This simply shows that the function v has a unique zero. That is 

, n , , ( b i + b 2)Y 7 P a 3*f> 

v(0) = a t +--- ( -. 

b 4 a ± ii 

From here, we observe that if a 4 + ^ bl+l>2 ^ Y > £ + we conclude that r(0) > 0 and thus, A 3 < 0 and this 

D4 CL'i [A 

implies that the steady state ^/ ai , 0, Vh 4 ^ ’ s asymptotically stable for all r > 0. If a 4 + 

we then conclude that v(0) < 0 and therefore A 3 > 0 which shows that the steady state ^/ ai , 0 J ^ is 
unstable for all r > 0. 


Notice that these roots changes as r varies. Now, we determine whether there exists a critical delay r c > 0 
interms of d t , d 2 , d 3 , k lt k 2 , k 3 so that Re(A ) > 0 for r > r c . This critical value of r is the point at which (3.10) 
transform from having a non positive to a non negative real parts. In other words, r c is the value of r such that 
Re (A) = 0, at which the transition from stability to instability occurs. By determining the conditions on the 
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parameters to ensure that the steady state ^/ ai , 0, ^/{j^j of the delay system (3.2) is still stable, we consider 
equation (3.10) as a complex variable mapping problem. 


For this to occur considering the steady state 



, there must be a boundary case, a critical value of 


r such that the characteristic has a purely imaginary root. To determine whether or not such a r exists, we 
reduce (3.10) to a polynomial problem and then seek for particular types of roots, thereby determining if a 
bifurcation can occur at such delay. We begin by looking for purely imaginary roots, and letting A(t) = 0(r) + 
iff (r), where 0 amd o are real, we have 0(0) < 0. By continuity of 0,0(r) < 0 for values of r such thatO < 
r < r c for some r c > 0. Therefore, the steady state remains stable for values of r. Substituting A(r) into 
equation (3.10) i.e., 


2A + diA 2 + d 2 X + d 3 + (A^A 2 + k 2 A + k 3 )e r ^ — 0 

we obtain 

(0 + id) 3 + di(0 + id) 2 + d 2 (0 + id) + d 3 

T (.k± (0 + iff) 2 + k 2 (0 + io - ) + k 3 )e~ r ^ +ia ^ = 0 

0 3 + 3<p 2 id — 3<pd 2 — id 3 + d^tp 2 + 2<pid — d 2 ) + d 2 (0 + id) + d 3 
+ [ki(0 2 + 2 <pid — d 2 ) + k 2 ((p + id) + k 3 ]e~ r ^ +ia ^ = 0 

0 3 + 3 (p 2 id — 3(pd 2 — id 3 + d 1 0 2 + 2d ± (pid — d x d 2 + d 2 <p + d 2 id + d 3 
+ [/c 1 0 2 + 2k 1 cpid — k ± d 2 + k 2 <p + k 2 id + k 3 ]e~ r ^ +lcr ^ = 0. 

We probably write the exponential in terms of trigonometric function, that is 

e -T(0+icr) _ e -T(t> e -Tia' 

But e~ Tl = cos Td — isin Td, 

=> e~ T< ^.e~ xia = e~ x< ^{cosTd — isinzd). 

Then we break the polynomial up into its real and imaginary parts, and write the exponential in terms of the 
trigonometric functions to get 

d 3 + d 2 (p + d^ 2 + 0 3 — d±d 2 — 3 <pd 2 + i(d 2 d + 2d 1 (pd + 3 (p 2 d — d 3 ) 

+e _T ^(cosT<T — isin rcr)[k 1 0 2 + 2/t!0iff — k x d 2 + k 2 (p + k 2 id + k 3 ] = 0 
d 3 + d 2 (p + d x 0 2 + 0 3 — d t d 2 — 3<pd 2 + i(d 2 d + 2d 1 cr + 3 <p 2 d — d 3 ) 

+e~ r( ^[cosTd (ki0 2 + 2 k ± (pid — k x d 2 + k 2 (p + k 2 id + k 3 )] 

— e~ r< ^[isin rcr(/c: 1 0 2 + 2k ± cpid — k x d 2 + k 2 (p + k 2 id + k 3 )] = 0 


d 3 + d 2 <p + d x 0 2 + 0 3 — d ± d 2 — 3<pd 2 + i(d 2 d + 2 d ± <pd + 3 <p 2 d — d 3 ) 
+e~ I( f > [(k 3 + k 2 (p + k x 0 2 — k-^d 2 ) cos Td + i{k 2 d + 2k t d) cos id] 
+e~ T( l > [(—k 1 d 2 + k 3 + k 2 (p + 0)(— isinrd) + i(k 2 d + 2k 1 cpd)(—isin Td)] = 0 

d 3 + d 2 (p + d!0 2 + 0 3 — d-^d 2 — 3(pd 2 + i(d 2 d + 2d 1 (pd + 3 (p 2 d — d 3 ) 
+e~ r< ^ [( k 2 d + 2k 1 (pd)sin Td + ( k 3 + k 2 (p + k 1 0 2 — /qcr 2 ) cos Td] 

+ie _T< ^ [(/tiff 2 — k 3 — k 2 (p — kiCp^sin Td + ( k 2 d + 2k ± (pd) cos tct] = 0. (3.11) 
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Assuming, we let 0(r c ) = 0 for some r c > 0 and 0(r) < 0 for 0 < r < t c , then the steady state 

may loss stability at r = r c or A(r c ) = io( t c ). In fact io is a root of equation (4.26) if and 

only if we take terms of o only and the constants. 

—ter 3 — d x o 2 + d 2 io + d 3 + (— k x o 2 + k 2 io + k 3 )(cos to — isin to) = 0 (3.12) 

In order for (3.12) to hold, both the real and imaginary parts must be zero. By equating real parts and imaginary 
parts of the right side of system (3.12) to zero, we have the equations; 

k 2 osin to + (/t 3 — k-^o 2 ) cos to — d ± o 2 + d 3 = 0 
k 2 io cos to — (k 3 — k-iO^isin to — io 3 + d 2 io = 0 



which gives; 


k 2 osin to + (fc 3 — k x o 2 ) cos to = d x o 2 — d 3 
k 2 io cos to — (k 3 — k ± o 2 )isin to = io 3 — d 2 io 


(3.13) 

(3.14) 


Squaring both sides of (4.29) and (4.30) one obtains; 

C k 2 osin.TO ) 2 + (k 3 — k^ 2 ) 2 (cos to) 2 = (^a 2 — d 3 ) 2 
=> k 2 o 2 sin 2 TO + {k 3 — 2k 1 k 3 o 2 + klo 4 )cos 2 TO = dfer 4 — 2d ± d 3 o 2 + d 3 


( k 2 ocos to) 2 + (k 3 — k 1 cr 2 ) 2 (—sinrer) 2 = (<j 3 — d 2 o) 


=> k 2 o 2 cos 2 To + (k 3 — 2k 1 k 3 o 2 + klo 4 )sin 2 TO = o 6 — 2d 2 o 4 + d 2 o 2 

Adding up the squares of (3.15) and (3.16) we have 

k 2 o 2 (sin 2 TO + cos 2 to) + [(k 2 — 2k 1 k 3 o 2 + klo 4 )(sin 2 TO + cos 2 to )] 

= d 2 o 4 — 2d 1 d 3 o 2 + d 3 + o 6 — 2 d 2 o 4 + d 2 o 2 

But sin 2 TO + cos 2 to = 1, then we have-, 

k 2 o 2 + k 2 — 2k 1 k 3 o 2 + k^o 4 = d\o 4 — 2d 1 d 3 o 2 + d 2 + o 6 — 2 d 2 o 4 + d\o 2 
=> U(o) = o 6 + (df — 2d 2 — k\)o 4 + (d| — 2d 1 d 3 + 2k t k 3 — k 2 )o 2 + d 2 — k 2 = 0 (3.17) 

After squaring (3.13) and (3.14) and then summed (3.15) and (3.16), we observe two things. The first 
observation is that the trigonometric term disappears from the polynomial and the delay, t, has been eliminated. 
The second observation tells us that the equation is an even polynomial. This is because when an even or odd 

function is squared, it will always result to an even polynomial, that is s(— o) 2 = (±s(er)) = s(cr) 2 
For further simplification of (4.33), we define a new variables by letting ; 


(3.15) 

(3.16) 


(, g = o 2 , h = di — 2 d 2 — k 2 , 1 = d 2 — 2 d x d 3 + 2k ± k 3 — k 2 , 

Then we have equation (4.33) written as 

N(g) = g 3 + hg 2 + Ig + m 1 = 0 


m 1 = d 2 — kl) G 
(3.18) 


where N is a polynomial. Remember that we are interested in o E R, therefore, if all the real part of N are non 
positive, it simply implies that there can be no simultaneous solution of (3.13) and (3.14). Alternatively, if there 
is a non negative real root g to N, there is a delay r corresponding to o which solves equations (3.13) and 
(314). 

Lemma 3: 

Assume that equation (3.18) has negative real roots, then all roots of equation (3.10) have negative real parts if 
r > 0. 


@ IJTSRD | Available Online @ www.ijtsrd.com | Volume-2 | Issue-5 | Jul-Aug 2018 


Page: 1686 



International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 

Proof: 

As far as equation (3.18) has negative real roots, any real number a is not a root of equation (3.17). This 
ensures that any real number a is not a root of equation (3.12). Hence for any real number a, the value io is not 
a root of equation (3.10) and this implies that there is no r c such that A(r c ) = io(x c ) is a root of equation (3.11). 
Because of proposition 2 , for a root of (3.10) corresponding to r = 0 , Re(A( 0 )) < 0 . Since i?e(A(r)) is a 
continuous function of x, we conclude that all roots of (3.10) have negative real parts. Since (3.18) is odd (of 
degree three), we are guaranteed a negative real root. The only way to have a simple positive real root in this 
case is to have two positive real roots. 

To present the conditions that ensure equation (3.18) has a positive real root or has no positive real roots, we 
take the Sturm chain of the polynomial (3.18) denoted by 

N 0 (g) = N(g); 
iVi(£) = N'(g). 

Therefore, N 0 (g) = N(g) = g 3 + hg 2 + Ig + 7 % = 0 (3.18) 

N ± (g) = N' (g) = 3g 2 + 2hg + l = 0 (3.19) 

The roots of (3.19) are 

_ -h + V/i 2 - 3 1 _ -h - V/i 2 - 31 

9i g ’ 9i ^ 1 

One of these is positive if h < 0 or h > 0 and l < 0, so either h or l must be negative. So we have; 

Lemma 4: 

A steady state with characteristic equation (3.18) is stable in the absence of delay and becomes unstable with 
increasing delay if and only if h, l and m are not all positive and 

1. d-L + k 4 > 0, c ?3 + k 3 > 0, (d 4 + k 4 )(d 2 + /C 2 ) — (c ?3 T k^) > 0 

2. Suppose that either (a) m < 0, or (b) m > 0, h 2 — 3/ > 0, h < 0 and Nig-t) < 0 then equation (3.18) has 
positive real roots. 

3. Suppose m > 0 and h 2 — 31 > 0, then equation (3.18) has no positive real roots. 

4. 

Proof: 

1 . d 1 = 2a 1 + b A +{b 1 + b 2 )^~, 

CLill 

d 2 = (% + ( b x + b 2 ) (% + h 4 ) + a 1 b 4 , 

k 2 = ~^ + ^ 7 ) Oi + b *)> 

d 3 = a 1 b 4 (a 4 + (b 4 + b 2 ) £), 
k 3 = -aA(( + 77> 

Suppose that the Routh-Hurwitz criterion holds. That is if d 4 + k 4 > 0, d 3 + k 3 > 0 and (d 4 + k 1 )(d 2 + 
k 2 ) — (d 3 + fe 3 ) > 0, then the characteristic equation (3.18) is stable in the absence of delay. 

(i) If (a) is true, that is if m < 0, then from equation (3.18) we have that A( 0 ) = m < 0. Also, we observe 
that 

lim N(g) = 00 , 


and this implies that the intermediate value theorem, guarantees that equation (3.18) must have a positive real 
root t 0 , with IV (t 0 ) = 0. 
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2. Now, suppose that (b) is true, that is if m > 0 and h 2 — 31 < 0 and h > 0, we find that g 1 is real and 
g 1 > 0. Since N( 0) = m > 0 and N(gf) < 0, also by the intermediate value theorem, N has a zero 
between the origin and g 1 . 

Notice that N'(g) = 0 has no real root since both zeros g 1 and g 2 of N' are not real if 3Z > h 2 . 

Therefore, 

h 2 

IV'(0) = l> — > 0 


In conclusion, we have that the coefficients; h, l and m of the quadratic polynomial N' are all positive real 
numbers and it simply means that N is increasing on the real numbers. Actually, because IV(0) = m > 0, we 
observe that when m > 0, N(g) does not vanish and thus, equation (3.18) does not have any positive real root. 
Also, lemma 4(ii) simply tells us that there is a such that a is a solution of the characteristic equation (3.10), 
thus indicating that real parts of all the eigenvalues of H are non positive for all delay r > 0. The summary of 
the above analysis is represented in the theorems below. 


Theorem 5: 

The uninfected steady state > 0, ^/b^j delay system (3.2) is asymptotically stable for all r > 0 if 

<■ + < a ± + ai ,a 3 , b 1 ,b 2 ,b 4 ,Y,^,f,^,g > 0, m > 0 and h 2 — 31 < 0. 


axil 


b 4 


Proof: 

By part (ii) of lemma 4, equation (3.18) has no positive roots. Lemma 3 ensures that all roots of (3.10) have 

negative real parts for r > 0. Therefore, the steady state 0,^/^ ^ of the delay system (3.2) is 

asymptotically stable for all r > 0. We will introduce the Hopf Bifurcation theorem (Hassard et al, 1981), as 
shall be discussed below. 

To establish Hopf bifurcation at r = t c , we need to show that 

— 0(r)^O. (3.20) 


Now to ensure that Hopf bifurcation occurs, we provide conditions on the parameters. Suppose the conditions 
in Lemma 4(ai) hold, then equation (3.17) has a positive root. We denote by gj,j E (1) or j E {1,2,3} 
depending on the number of positive roots that equation (3.18) has. Observe that equation (3.17) has up to six 
roots i.e., ±Jg] for j = 1,2,3. 


Notice that if the solution of equation (3.9) exists, it is among these +yfg] for j = 1,2,3 and that if A = io is a 
root of equation (3.10), —ia will also be a root. 


To equations (3.13) and (3.14) we substitute <t ; - and then solve for r, which will then be of the form 


(n) L ■ - 

t) = — sin 


kitf + (d 4 k 2 — k 3 — d 2 k 4 )af + (d 1 /c 3 — d 3 k 2 )af 

L 

2n(n — 1) 

k 2 a 2 + (k 3 - k 4 a 2 ) 2 

1 

. °j . 


where j = 0,1,2. 


if we multiply (3.13) by k 2 a c and (3.14) by k 3 — k t a 2 and also add the resulting product, we obtain 

k 2 o 2 sin o c t c + k 2 a c (k 3 — k 4 a 2 ) cos <t c t c = k^^d ^ 2 — d 3 ) 

-k 2 a c (k 3 - k 4 a 2 ) cosa c r c + ( k 3 - k ± a 2 ) 2 sina- c r c = ~(k 3 - k ± a 2 )(a 2 - d 2 o c ) 
=> \kW + (k 3 - k 4 a 2 ) 2 ] sin<r c T c = k^dd-^a 2 - d 3 ) - (k 3 - k 4 a 2 )(a 2 - d 2 o c ) 
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_ k 2 a c {d x a 2 - d 3 ) - (fc 3 - k t a 2 )(a 3 - d 2 a c ) 
sm (T c T c - [fc^ + fe-Vc 2 ) 2 ] 


1 ■ -1 
x c = — sin 1 

kiOf + (d 1 /c 2 — k 3 — d 2 ki)a 2 + ( d 2 k 3 — d 3 k 2 )af 

+ 

2n(n — 1) 

k 2 a 2 + (k 3 - k t a 2 ) 2 

a J . 

x\ n = — sin -1 

J Gj 

'k ± af + (d ± k 2 — k 3 — d 2 k-P)af + ( d 2 k 3 - d 3 k 2 )af 

+ 

2n(n — 1) 

k\aj + (fc 3 - k x a 2 ) 2 

. Gj 


Now let r c be the smallest of such x for which p(x c ) = 0. 


r =r (nc) 
Therefore, c 


= min{r| n) > 0,0 < j < 2,n > l}| ^ ^ 

^ = ° jc J 


Theorem 6: 

Let the critical time lag x c and a c be defined as (3.21) for the time lag x and suppose that 3 ap + 2her,? + /a,? =£ 
0 then the system of delay differential equation (3.2) exhibits the Hopf bifurcation at the steady 


state! ^/ ai ,0//,|. 


Proof: 

To establish Hopf bifurcation, we show that 

dp(x) 


(3.22) 


and this guarantees that Hopf bifurcation occurs. For (3.21) to hold, both the real and imaginary parts of (3.11) 
must be zero. So we get the pair of equations 

e~ T [(.k 2 a + 2/r 1 0cr) sin (pa + ( k 3 + k 2 (p — k t p 2 — /qcr 2 ) cos pa] 

+d 3 + d 2 p + d t p 2 + p 3 — d±a 2 — 3 pa 2 = 0 
ie~ r [(/cya 2 — k 3 — k 2 p — kip 2 ) sin pa + ( k 2 a + 2/c 1 0a) cos pa] 

+i(d 2 a + 2 d ± pa + 3 p 2 a — a 3 ) = 0 


(3.23) 

(3.24) 


which finally gives 


e~ r $ [(/c 2 a + 2/c 1 0cr) sin pa + ( k 3 + k 2 p — kip 2 — k-^a 2 ) cos pa] 
= (c?! + 30)cr 2 — d 3 — d 2 p — dxp 2 — 0 3 

e _T ^ [(/cycr 2 — k 3 — k 2 p — kip 2 ) sin pa + ( k 2 a + 2k ± pa) cos pa] 
= a 3 — d 2 a — 2dipa — 3 p 2 a 


(3.25) 

(3.26) 


From equation (3.25), we have that e~ T( ^ = e° since xp = 0. Then e° = 1, and selecting a from the equation 
we obtain 

k 2 a sin pa + k 3 cos pa — k-^a 2 cos pa = g^ct 2 — d 3 (3.27) 


We at this moment differentiate equation (3.27) with respect to x after which we evaluate it at r = x c for which 
0 (t c ) = 0 and a(x c ) = a c and get 
da 

— | T=T = k 2 sin a c x c + x c k 2 a c cos a c x c — k 3 x c sin a c x c — 2k t a c cos a c x c + x c kia c sin a c x c 


dx 


— 2d 1 cr c 

2d 1 cr c + (2 kia c — x c k 2 a c ) cos a c x c + ( k 3 x c — x c kia c 2 — k 2 ) sin a c x c = A x 
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Select 0 from (3.25), differentiate with respect to r and evaluating at r = r c for which 0(r c ) = 0 and <t(t c ) = 
a c , we then obtain 

2/c 1 0cr sin 0(7 + k 2 (p cos 0<7 + fci0 2 cos 0<7 = 3 (pa 2 — d 2 (p — d t (p 2 — 0 3 (3.28) 


d0 

dr h 


But 0(T c ) = 0 


2/cycr sin <7 c t c +2r c 0 /qcr,- cos <7 c t c + /r 2 cos <j c t c — x c k 2 (p sin <7 c r c + 2 k 4 (p cos <7 c r c 
— k 1 (p 2 x c sin <7 c t c = 3<7 c 2 — d 2 — 2d 1 0 — 30 2 

=> 3<7 c 2 — d 2 — 2k x a sin <7 c r c — /c 2 cos <7 c r c = A-l 


Also from equation (3.26) we have e T< ^ = e° = 1. Selecting <7 from (3.26) and differentiating with respect to x 
and evaluating at r = r c for which 0(r c ) = 0 and a(x c ) = a c we get 

k r a 2 sin ax — k 3 sin ax + k 2 a cos ax = a 3 — d 2 a (3.29) 


Then, 


da 

dx It=Tc 


= 2k 1 a c sin a c x c + k-^Xf-a 2 cos a c x c — x c k 3 cos a c x c — k 2 x c a c sin a c x c + k 2 cos a c x c 

= 3a c 2 — d 2 

d 2 — 3a c 2 + (k 2 + kxXca 2 — x c k 3 ) cos a c x c + (2/c 1 <t c — k 2 x c a c ) sin a c x c = A 2 


Again, selecting 0 from (3.26), then and differentiating with respect to r and evaluating at r = x c for which 
0(t c ) = 0 and a(x c ) = a c we have 

—/c 2 0 cos ax — /Ci0 2 sin ax + 2/c 1 0a cos ax = — 2d 1 0<7 — 3 (p 3 a (3.30) 

Then, 

d(p 2 

— | T=Tc = ~k 2 sin a c x c - k 2 x c (p c cos a c x c - 2k 1 (p c sin a c x c - k ± x c (p c sin a c x c 

3~2k 4 a c cos a c x c — 2ki<f> c a c sin a c x c = — 2d-±a c — 6<p c a c 
=> —k 2 sin a c x c + 2k 1 a c cos a c x c = —2d 1 cr c 
=> 2d 1 <7 c — k 2 sin <7 c r c + 2k 1 a c cos a c x c = A 2 

At this point, we solve for the difference of A 1 and A 2 and the sum of A 2 and A 1 , that is 

A i ^ \t=t c ~ A 2^ \t=t c and A 2 ^| T=Tc + Ai ^ \ T=Tc respectively, then we obtain 

da dtp 

i4l d7 lT=Tc_i42 d7 lT=Tc 

=> [2d 1 a c + (2 k^c — x c k 2 a c ) cos a c x c + ( k 3 x c — x^a 2 — k 2 ) sin a c x c ] 

—[2d 1 a c — k 2 sin <7 c r c + 2k 1 a c cos a c x c ] 

=> —x c k 2 a c cos a c x c + ( k 3 x c — x^a 2 ) sin a c x c 


Sincer c = a c , we have 


=> k 2 a c 2 cos a c x c + (k-ta 2 — k 3 a c ) sin a c x c 


Therefore, 

Also 


A 3 = k 2 a c 2 and A 4 = k^ 2 — k 3 a c 


da 

dx 


dd> 

+ A i~r:U 


A 2 ^ lr=r c ' dT It=t c 


[d 2 — 3cr c 2 + (k 2 + k 4 x c a c 2 — x c k 3 ) cos a c x c + (2/c 1 cr c — k 2 x c a c ) sin <7 c t c ] 
+ [3<7 c 2 — d 2 — 2k 1 asin a c x c — k 2 cos <7 c t c ] 
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Again remember that x c = a c , and we then obtain 

—k 2 a c 2 sin a c r c + (/cycr ,- 3 — k 3 a c ) cos <t c t c 
=> k 2 a 2 sin a c x c - (k 4 a c 3 - k 3 a c ) cos a c x c 


Therefore, 


A 3 = k 2 a c 2 and — A 4 = ( k x o 3 — k 3 o c ) 


The result of the differentiations, the difference and sum when put together, yields 
da d(p 

A 1 — \ T=Tc - A 2 — \ T=Tc = A 3 cosa c x c + A 4 sin a c x c (3.31) 

da d<p 

A 2 \r=r c + \r=r c = A 3 sin a c t c - A 4 cos a c T c (3.32) 

To solve equations (3.31) and (3.32) simultaneously, we multiply equation (3.31) by A 2 and equation (3.32) by 
A t , then subtract the product of (3.31) from (3.32) and this eliminates ^ .Therefore, we have 


da 


l2 d <t> 


A X A 2 It ' T=Tc + Al ~dx ' T=Tc = AlAs Sin ffcTc ~ A±A4 C ° S acTc 


da 


2 dtp 


A 4 A 2 —— |x=r — A 2 — | t=t = A 2 A 3 cos a c T c + A 2 A 4 sin a c r c 


dr 


d(b 

(A 2 + A\) — | T=Tc = (A 4 A 3 - A 2 A 4 ) sin a c x c - ( A t A 4 + A 2 A 3 ) cos a c x c 

(3.33) 


dx 


d4 

dx It=Tc 


dx IT=Tc 

( A t A 3 — A 2 A 4 ) sin a c x c — (A 1 A 4 + A 2 A 3 ) cos a c x c 


A 2 + A 2 


Let us solve equations (3.13) and (3.14) simultaneously. 

k 2 a c sin a c x c + ( k 3 - k 4 a 2 ) cos a c x c = d 4 a 2 - d 3 

k 2 a c cos a c x c — ( k 3 — k ± a 2 ) sin a c x c = a 3 — d 2 a c 


(3.13) 

(3.14) 


If we multiply (3.13) by k 3 — k t a 2 and (3.14) by k 2 a c and add the resulting products, we obtain 
k 2 a c (k 3 - k 4 a 2 ) sin a c x c + ( k 3 - k 4 a 2 ) 2 cos a c x c = (d 4 a 2 - d 3 )(k 3 - k 4 a 2 ) 

k 2 a 2 cos a c x c - k 2 a c (k 3 - k 4 a 2 ) sin a c x c = k 2 a c (a 3 - d 2 a c ) 

=> kja 2 cos a c x c + (k 3 - k 4 a 2 ) 2 cos a c x c = (d 4 a 2 - d 3 )(k 3 - k 4 a 2 ) + k 2 a c (a 3 - d 2 a c ) 


\k\a 2 + (k 3 - k 4 a 2 ) 2 ] cos a c x c = (d 4 a 2 - d 3 )(k 3 - k 4 a 2 ) + k 2 a c {a 3 - d 2 a c ) 


kiVc&l ~ g c ) + ( 'Aoc ~ d 3 )(fc 3 ~ feig c 2 ) 

\.k 2 a 2 + (k 3 - k 4 alY] 


(3.34) 


On the other hand, if we multiply (3.13) by k 2 a c and (3.14) by k 3 — k 4 a 2 and also add the resulting product, 
we obtain 

k 2 a 2 s\na c x c + k 2 a c (k 3 - k 4 a 2 ) cosa c x c = k^dd^ 2 - d 3 ) 


-k 2 a c (k 3 - k 4 a 2 ) cosa c x c + (k 3 - k 4 a 2 ) 2 sina c x c = ~(k 3 - k ± a 2 )(a 3 - d 2 a c ) 


=> [kW + (k 3 ~ k ± a 2 ) 2 ] sin<j c r c = k 2 a c {d 4 a 2 - d 3 ) - (k 3 - k 4 a 2 )(a 3 - d 2 a c ) 
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k 2 (J c (d 1 a 2 - d 3 ) - (fc 3 - /ti<7 c 2 )((X c 3 - d 2 o c ) 


[/c|a c 2 + (k 3 - /qa 2 ) 2 ] 


(3.35) 


Substituting (3.34) and (3.35) into equation (3.33) it yields (3.36) below 
deb 3(Jr + 2 h.Or + la} 

-j- \ T - T =— - 5 - ^0 

dr c A\+A\ 


(3.36) 


Therefore we conclude that Hopf bifurcation occurs when r passes through the critical value which is assumed 


to be r c . 


IV. ANALYSIS OF RESULT 

Certainly, our interest is in determining if the steady 
state of the model system is stable or not. When the 
slightest disturbance completely alters the behavior of 
the system, them it is impossible to predict the long 
term evolution of the system under consideration. We 
therefore showed in this research work with the 
parameter value by Long et al (2008) that critical 
delay beyond which the steady state would lose its 
stability through oscillation does not exist. 

Biologically, proposition (4i) states that if a 1 + 
( 'b 1 +b 2 )Y t [ ien q ie steady state (—, 0, —) 

b 4 Vax b 4 y 

of system (4.18) is asymptotically stable and this 
implies that there will be viral clearance because, 
CTL will both cytolitically and non cytolitically take 
care of infected hepatocytes before they can produce 
sufficient amounts of viral progeny which finally 
results in the absence of the virus. This is shown in 
figure (la - lc) and (2a - 2c). Figure Id focused on 
the behavior of the density of hepatitis B virus 
infected hepatocytes and CTL when the delay 
parameter r = 1 and show the interaction between 
virus and immune response. 






Figure (la - Id): Numerical simulation of the delay 
model system (18) where parameter values, a l = 
0.002, y = IF — 3, a 3 = 1500, a 4 = 0.3, b 1 — 
500, b 2 = 1500 ,x/j = 1 ,b 4 = 0.2, £ = 0.04,^ = 


0.58, [3 = 10 (Long et al,200S)and r 

= 1 (Eikenberry et al, 2009). 




Figure (2a - 2c): Numerical simulation of the delay 
model system (18) where parameter values, a 1 = 
0.002,y = IE — 3, a 3 = 1500, a 4 = 0.3, b 4 = 
500 ,b 2 = 1500,1/; = 1 ,b 4 = 0.2, ( = 0.04,^ = 
0.58, (3 = 10 and r = 1 with time ranging from 0 to 

2.5. 



Figure (3a - 3c): Numerical simulation of the delay 
model system (18) where parameter values, a t = 
0.002,y = IE — 3, a 3 = 1500, a 4 = 0.3 ,b 4 = 
500, b 2 = 1500 ,ip = 1 ,b 4 = 0.2, ( = 0.04,^ = 
0.58, (3 = 10 and x = 1 with time ranging from 0 to 

50 
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Figure (la - Id), (2a - 2c) and (3a - 3c) shows the 
delay effect on equation (3.2) with delay value r = 1. 
The infected steady state of proposition (li) of delay 
system (3.2) remains stable if the conditions in 
theorem 5 which states that the uninfected steady state 

of the delay system (3.2) is 

asymptotically stable for all r > 0 if ( + < a 4 + 

(fcl ^ 2)r - a 3 > K b 2 , b 4 , y, /?, ip, p > 0, m > 

0 and h 2 — 31 < 0 

are satisfied. Then, we can conclude that 
incorporating a discrete intracellular time delay, into 
system (3.1) does not change the stability of the 
steady state of proposition (li). In this case, even if 
the immune response is maintained and virus persists, 
the virus may persist without causing much harm to 
the liver since the immune system is strongly 
responsive in eliminating the virus from within the 
hepatocytes without killing the hepatocytes. 

If theorem 8 which says that suppose that 3 of + 
2h<Jc + la 2 =£ 0 then the system of delay differential 
equation (4.18) exhibits the Hopf bifurcation at the 

steady state (~’0, is satisfied, we have that the 

stability of the infected steady state of proposition (li) 
of delay model equation (3.2) depends on the delay 
value and the delay can induce oscillation in the 
number of the uninfected hepatocytes, infected 
hepatocytes and CTL. If the individuals’ ALT level is 
high, then, the virus may have escaped the control of 
the CTL and HBV infection is reactivated. However, 
if the immune control succeeds, the virus population 
will utterly be destroyed. Ciupe et al (2007), 
demonstrated that consideration of delayed cytotoxic 
and non cytotoxic immune reactions and the presence 
of cells refactory to infection was necessary to 
properly understand the dynamics of HBV acute 
infection and progression to chronic disease. 

V. SUMMARY AND CONCLUSION 

In summary, we have that clearance of acute HBV 
infection occurs subsequent to the appearance of the 
Cytotoxic T Lymphocyte (CTL) cells. The clearance 
may be cytolytic or noncytolytic. These events are 
accompanied by a significant increase in apoptosis 
and regeneration of hepatocytes as indicated in figure 
(2a - 2c) and (3a - 3c). the liver has a mark ability to 


regenerate. Resection of half of the human liver is 
followed within 2-4 weeks by full restoration of the 
liver structure, size and function. The influx of the 
cytotoxic T lymphocytes will delay the accumulation 
of virus free hepatocyte. More so, chronically infected 
individuals can exhibit to some level T - cell 
accumulation, cytokine expression, and apoptosis. 
Our model predicted replacement of infected 
hepatocytes with regenerated cells which by known 
mechanism remain protected from reinfection in the 
individual after immune response. 
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